Hole dynamics in noble metals 
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We present a detailed analysis of hole dynamics in noble metals (Cu and Au), by means of first- 
principles many-body calculations. While holes in a free-electron gas are known to live shorter than 
electrons with the same excitation energy, our results indicate that d-holes in noble metals exhibit 
longer inelastic lifetimes than excited sp-electrons, in agreement with experiment. The density of 
states available for d-hole decay is larger than that for the decay of excited electrons; however, the 
small overlap between d- and sp-states below the Fermi level increases the d-hole lifetime. The 
impact of d-hole dynamics on electron-hole correlation effects, which are of relevance in the analysis 
of time-resolved two-photon photoemission experiments, is also addressed. 



Detailed and quantitative understanding of hot- 
electron and photo-hole dynamics is the prerequisite for 
any tailoring of technological processes in solid-state 
physics and surface chemistry, which are governed by 
charge transfer and electronic excitations In par- 

ticular, the quasiparicle lifetime, i.e., the time a quasipar- 
ticle propagates without loosing its phase memory, rep- 
resents a basic quantity in the analysis of these processes. 
Experimentally, angle-resolved photoemission (PE) spec- 
troscopy provides direct information on hole lifetimes 
A new path for the study of both electron and 
hole (and spin) dynamics in the time domain was opened 
by the advent of the time-resolved two-photon photoe- 
mission (TR-TPPE) technique ^j^. In these "pump- 
probe" experiments, the emitted photoelectron is mea- 
sured both with energy and momentum resolution; elec- 
tron and hole lifetimes are then measured at well-defined 
k-space points, by combining the band mapping capabil- 
ities of photoemission with the time resolution of nonlin- 
ear optical spectroscopy [ pT| . 

The theoretical framework to investigate the quasipar- 
ticle lifetime has been based for many years on the free- 
electron gas (FEG) model of Fermi liquids ||l3L char- 
acterized by the electron-density parameter rg |13[| . In 
this simple model and for either electrons or holes with 
energy E very near the Fermi level {E ^ Ep), the in- 
elastic lifetime is found to be, in the high-density limit 
(r, « 1), t{E) = 263r~^^^ (E- EF)-^fs, where E 
and Ep are expressed in eV Several other free- 

electron calculations of electron-electron scattering rates 
have also been carried out, for electron/hole energies that 
deviate from the Fermi level, within the random-phase 
approximation (RPA) ||l^,|l^ and with inclusion of ex- 



change and correlation effects |l7| , |l8| . Nevertheless, de- 
tailed TR-TPPE experiments have reported large devia- 
tions of measured hot-electron lifetimes from those pre- 
dicted within the FEG model ||l|-|||l. Moreover, while 
holes {E < Ep) in a FEG are known to live shorter than 
their electron counterparts {E > Ep) [for the same exci- 
tation energy, the momentum of the hole is smaller than 
the electron momentum, thus yielding a larger number 
of available transitions], recent TR-2PPE measurements 
[p6| have demonstrated clearly that the d-hole lifetime 
at the top of the Cu d-bands is considerably slower than 
that of excited electrons with the same excitation en- 
ergy. Subsequently, recent PE measurements have 
arrived at the same conclusion for all three noble metals 
Cu, Ag, and Au. First-principles calculations of elec- 
tron lifetimes in a variety of metals have been carried 
out only very recently |^,^, by using the GW approx- 
imation of many-body theory |Q. These calculations 
pl| show that band-structure effects play a key role in 
the quasiparticle-decay mechanism. Furthermore, hot- 
electron lifetimes are found to strongly depend on the 
momentum of the quasiparticle, especially in the case of 
metals with a band structure having van Hove-like sin- 
gularities in the vicinity of the Fermi level . 

During the electron-pump process induced by the first 
ultrashort laser in TR-TPPE experiments, a hole is cre- 
ated in the Fermi sea. This hole can be filled with an 
electron from below the Fermi level, via Auger decay, 
higher-order excitations, or electron-phonon interactions. 
The understanding of the hole decay is fundamental for 
the description of the ulterior dynamics of the excited 
electron (probed by the second laser). Furthermore, the 
dynamics of the hole is interesting by itself, as it is a 
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fingerprint of differences between electron and hole wave 
functions. These differences yield distinct behaviours of 
electron and hole lifetimes, which depend on both the 
energy and the momentum of the quasiparticle. 

In the present letter, we address the first stage of a TR- 
2PPE experiment by investigating the hole-quasiparticle 
dynamics, and report first-principles many-body calcu- 
lations of the decay rate of d-like holes in Cu and Au. 
Details of the formalism we use for the evaluation of de- 
cay rates are described in Ref. |2|]. The basic equation 
for the damping rate of a hole in the state 0n,k(r) with 
energy £„.k < Ep is (we use atomic units throughout. 
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where E(r, r';e„^k) represents the quasiparticle self- 
energy, which we compute within the GW approximation 
of many-body theory ||30(| . For the evaluation of both the 
initial state of Eq. (Ill) and all wave functions entering 
the quasiparticle self-energy, we first expand the one- 
electron Bloch states in a plane- wave basis (PW), and 
then solve self-consistently the Kohn-Sham equation |3|] 
of density- functional theory (DPT). Though all-electron 
schemes, such as the full-potential linearized augmented 
plane- wave (LAPW) method , are expected to be bet- 
ter suited for the description of d-bands, the PW pseu- 
dopotential approach has already been succesfuUy incor- 
porated in the description of the dynamical response of 
copper |3^. Both PW and LAPW calculations produce 
almost identical results, showing the correct overall band 
structure. Nevertheless, they both predict the entire d- 
band manifold to be ^ 0.5 eV higher than shown by pho- 
toemission experiments [ ^HIJb^ . 

The quasiparticle decay rate in periodic crystals de- 
pends on both the wave vector k and the band index 
n of the initial Bloch state. Nevertheless, we also de- 
fine T^^{E), as the average of T~^(k, n) over all wave 
vectors and bands lying with the same energy E in the 
irreducible wedge of the BZ. Our full band-structure cal- 
culations of the average lifetime t{E) of holes {E < Ep) 
in Cu and Au are presented by open circles in Figs, la 
and lb, respectively ||3^,Q. For comparison, full band- 
structure calculations of the average lifetime of electrons 
{E > Ep) are exhibited by solid circles, and PEG cal- 
culations [with rs — 2.67 for Cu and = 3.01 for Au] 
are represented by solid {E > Ep) and dotted {E < Ep) 
lines. These results indicate that holes in noble metals 
exhibit considerably longer lifetimes than electrons with 
the same excitation energy, which is in agreement with 
accurate TR-2PPE and PE measurements of d-hole life- 
times reported in Refs. |2^] and ||2^, respectively. More- 
over, at energies near the top of the d-bands the hole 
lifetime strongly deviates from the t{E) oc {E — Ep)'"^ 



quadratic dependence predicted within the FEG model, 
as first pointed out by Petek et al 

Both electrons and holes exhibit lifetimes that are well 
over those predicted within the FEG model of the solid, 
due to a major contribution from occupied d states par- 
ticipating in the screening of electron-electron interac- 
tions, and differences between electron and hole lifetimes 
stem from the intrinsic properties of Bloch states above 
and below the Fermi level. Whereas Bloch states above 
the Fermi level have mainly sp-like character, states just 
below the Fermi level have a small but significant d- 
component, which increases when the d-band threshold 
is reached. This difference accounts for the larger hole 
lifetime for energies very near the Fermi level, as shown 
in Fig. 1. At the top of the c?-bands [1.5 and 1.7 eV be- 
low the Fermi level in Cu and Au, respectively], the low 
overlap between d- and sp-states below the Fermi level 
yields a dramatic increase in the hole lifetime, especially 
in the case of Cu, which cannot be explained with use of 
the FEG model. As the hole energy decreases {\E — Ep\ 
increases), both the increased phase space for the hole 
to decay and the larger overlap for hole-hole scattering 
within the d-bands lead to a rapid decrease of the hole 



lifetime [26 1. 

This physical scenario for the hole dynamics in noble 
metals is drawn in Fig. 2, where the ratio between hole 
and electron lifetimes in Cu (solid circles) and Au (open 
circles) is plotted as a function of the hole/electron en- 
ergy with respect to the Fermi level. Though calculated 
lifetimes of both electrons and holes in Au are much 
longer than in Cu, due to a larger participation of oc- 
cupied d states in the dynamical screening of the former 
[ p9[ , different ratios between hole and electron lifetimes 
are only observed at energies near the top of the d bands. 
As 5c?-bands in Au are more free-electron-like than Sd- 
bands in Cu, the larger overlap between d and sp-states 
below the Fermi level in Au yields an increase in the hole 
lifetime at \E — Ep\ ~ 1.7 eV in this material that is 
not as dramatic as in the case of Cu. This different be- 
haviour of hole lifetimes in Cu and Au emphasizes the 
crucial role that the overlap of the d-hole with unoccu- 
pied sp- and d-hole states play in the hole decay mecha- 
nism. As \E — Ep\ increases, the ratio between hole and 
electron lifetimes resembles that predicted by the FEG 
model of the solid, which is plotted in the inset of Fig. 2 
for r,, = 2.67. 

d-bands in the noble metals open in the vicinity of 
the high-symmetry X-point. Hence, we have calculated 
lifetimes of holes with the wave vector along the TX di- 
rection in Cu. At the top of the d-bands, at the ^5 
point with E — Ep = — 1.5eV and the X2 point with 
E — Ep — —1.7, we obtain d-hole lifetimes of 99 and 90 fs, 
respectively, much larger than lifetimes of electrons with 
the same energy, showing the role that the low overlap of 
d-holes with unoccupied sp-hole states play in the deter- 
mination of the hole lifetime. Though these calculated 
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values of d-hole lifetimes cannot be directly compared 
with experiment, since the measured d-band threshold in 
Cu is located at ~ 2 eV below the Fermi level, our calcu- 
lations are in agreement with the experimental observa- 
tion that the narrowest linewidths correspond to the top 
of the c?-bands 1 26 2^ . In the case of holes with the wave 
vector along the TL direction in Cu, the hump of Fig. la 
at E — Ep ^ — 1.5eV is absent and the calculated life- 
times decrease with a quadratic [r cx [E — Ep)^"^] energy 
scaling. Instead, the hump of Fig. la is mainly origi- 
nated in the contribution to the average lifetime from k 
vectors along the TX direction, where the opening of the 
d-bands occurs, thereby showing a distinct behaviour of 
hole lifetimes along various symmetry directions ||40|] . 

As d-holes in noble metals are found to live longer than 
sp-electrons with the same excitation energy, the major- 
ity of electrons excited by the first probe pulse in a TR- 
2PPE experiment feels the field created by the existing 
hole, thereby altering the electron dynamics. Hence, the 
TR-2PPE hot-electron lifetime measurements include, in 
a complex way, contributions from the joint electron-hole 
(exciton) dynamics |^l[]. In the case of semiconductors, 
this excitonic renormalization is known to strongly mod- 
ify the single-particle optical absorption profile, so it cer- 
tainly needs to be included in electronic-screening cal- 
culations 1 ^^ . In principle, excitonic normalization does 
not play such an important role in the case of metals, be- 
cause of the large dynamical screening in these materials. 
However, there is a femtosecond time scale for the build- 
ing up of the screening. During this ultrashort time, the 
electron-hole interaction is not fully screened and might 
modify the lifetime of the excited electron, as measured 
by the second probe pulse. Hence, there is still much to 
be done to understand the time-dependent building up 
of the screening in real metals, as well as the many-body 
electron and hole correlation effects in the quasiparticle 
dynamics. 

In summary, we have presented a detailed theoreti- 
cal investigation of inelastic lifetimes of holes in Cu and 
Au. We have reached the important conclusion that d- 
holes in these materials exhibit longer lifetimes than sp- 
electrons with the same excitation energy, in agreement 
with experiment. While a major contribution from oc- 
cupied c?-states participating in the screening of electron- 
electron interactions yields both electron and hole life- 
times that are much longer than those of electrons and 
holes in a FEG, the small overlap between d- and sp- 
states below the Fermi level is responsible for the life- 
time of d-holes being longer than that of sp-electrons, 
especially at the top of the d-band. As the d-hole energy 
decreases {\E — Ef \ increases), both the increased density 
of states available for the hole-decay mechanism and the 
larger overlap for hole-hole scattering within the d-bands 
yield a rapid decrease of the hole lifetime. Although our 
full band-structure calculations predict longer lifetimes 
than measured in the experiment p3| , they are in agree- 



ment with the experimental observation that the narrow- 
est linewidths correspond to the top of the d-bands. Also, 
our results highlight new effects related to the electron- 
hole interaction occurring during the time delay between 
the two laser pulses in TR-2PPE spectroscopy, which 
may be of great importance in the interpretation of these 
experiments. The present results are general, they can 
be applied to the study of other d-metals, and can also 
be extended to the investigation of spin dynamics. 
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FIG. 1. Electron and hole lifetimes in (a) Cu and (b) An. 
Solid and open circles represent our full ab initio calculation 
of r{E) for electrons {E > Ef) and holes [E < Ep), respec- 
tively, as obtained after averaging r(k,n)-^ of Eq. (§) over 
wave vectors and the band structure for each k. The solid and 
dotted lines represent the corresponding lifetime of electrons 
(solid line) and holes (dotted line) in a FEG with = 2.67 
for Cu and — 3.01 for An. 

FIG. 2. Ratio between hole and electron lifetimes in Cu 
(solid circles) and Au (open circles). The inset exhibits the 
corresponding ratio between hole and electron lifetimes in a 
FEG with rs = 2.67. 
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Figure 1b 
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